bind
  lambda 10*(10^4)
  mu 1.25*10^6
end

markov M2(n,k)
loop i,n,1,-1
  * 0 to 1,  1 to 2 , ..., 
  $(i) $(i-1)    (i)*lambda
   $(i-1) $(i)    (n-i+1)*mu
end
reward default 1
loop i,0,k-1
  * declare states 0 to k-1 as down states  
  $(i)   0
end
end
* Initial probability
   $(n)  1
end


func A(n,k) exrss(M2;n,k)
func U(n,k) 1-A(n,k)
func MTTSF(n,k) A(n,k)/(k*lambda*prob(M2,$(k);n,k))
func MTTSR(n,k) U(n,k)/((n-k+1)*mu*prob(M2,$(k-1);n,k))
func MTTSR1(n,k) U(n,k)*MTTSF(n,k)/A(n,k) 


* vary the value k
bind n 4
loop k,1,3
  expr A(n,k), U(n,k), MTTSF(n,k), MTTSR(n,k), MTTSR1(n,k) 
end 

end




